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ABSTRACT 

Beginning in 1999 January, the bright, strongly variable Narrow-Line Seyfert 1 
(NLS1) galaxy Akn 564 has been observed by RXTE once every ~4.3 days. It was 
also monitored every ~3.2 hr throughout 2000 July. These evenly-sampled observations 
have allowed the first quantitative comparison of long and short time-scale X-ray vari- 
ability in an NLS1 and the derivation of an X-ray Power Density Spectrum (PDS). The 
variability amplitude in the short time-scale light curve is very similar to that in the 
long time-scale light curve, in marked contrast to the stronger variability on longer time- 
scales which is characteristic of "normal" broad-line Seyfert Is (BLSls). Furthermore, 
the Akn 564 PDS power law cuts off at a frequency of 8.7 x 10 -7 Hz corresponding to 
a timescale of ~13 d, significantly shorter than that seen in the PDS of NGC 3516, a 
BLS1 of comparable luminosity. 

This result is consistent with NLSls showing faster (as opposed to larger amplitude) 
variations than BLSls, providing further evidence that NLSls harbour lower mass black 
holes than BLSls of similar luminosity, accreting at a correspondingly higher relative 
rate. 

Subject headings: galaxies: active — galaxies: individual (Akn 564) — galaxies: Seyfert 
— X-rays: galaxies 

1. Introduction 

A relatively recent development in the taxonomy of active galactic nuclei (AGN) has been 
the emergence of Narrow-Line Seyfert 1 (NLS1) galaxies as an important sub-class. NLSls were 
originally identified on the basis of their optical properties: exhibiting narrower permitted lines (H/3 
FWHM < 2000 km/s) than "normal" Broad-Line Seyfert Is (BLSls) and weaker [O III]/H/3 < 3 
than Seyfert 2s (Osterbrock & Pogge 1985; Goodrich 1989) . ROSAT data showed these properties 
to be strongly correlated with X-ray spectral slope, in the sense that NLSls tend to have steep soft 
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X-ray spectra (Boiler et al. 1996). In addition, many NLSls are strongly variable on time-scales of 
hours or less, with some NLSls showing giant (factor of ~100) X-ray flares over time-scales of days 
(e.g., Brandt et al. 1999). On these short time-scales the variability levels in NLSls are typically 
a factor of ~3-10 times larger than those seen in BLSls of similar luminosity (Turner et al. 1999; 
Leighly 1999). 

A number of models have been proposed to explain the distinctive properties observed in 
NLSls. The most widely accepted explanation is that NLSls are accreting near the Eddington limit 
(Pounds et al. 1995, Laor et al. 1997), while BLSls accrete at a lower rate. On this interpretation 
the high accretion rate is directly responsible for the strong soft X-ray emission, by enhanced 
thermal radiation from the accretion disc, while the increased soft photon flux cools the hard X-ray 
source leading to the steeper power law frequently seen at higher X-ray energies in NLSls (Brandt 
et al. 1997). It is the prospect that comparative study of NLSls and BLSls will shed light on the 
accretion rate, one of the fundamental parameters of an AGN, that makes the emergence of the 
new sub-class of AGN of such interest and importance. 

X-ray variability has been recognised as a powerful probe of the central regions of AGN since 
the EXOSAT "long- looks" first showed (McHardy 1988) that rapid variability was common in 
Seyfert galaxies, lending early support to the now-standard black hole/accretion disc paradigm. 
The EXOSAT data were found to be well-described by a fluctuation power density spectrum (PDS) 
rising smoothly to lower frequencies as a power law, f~ a where a = 1 — 2, to a limit imposed by the 
maximum duration of the ~ 2-day observations (e.g., Lawrence & Papadakis 1993). Attempts to 
constrain a flattening or "cut-off' in the PDS (as required to avoid the variability power diverging) 
were made by combining EXOSAT data with earlier X-ray satellite observations of NGC 5506 
(McHardy 1988) and NGC 4151 (Papadakis and McHardy 1995). These analyses yielded evidence 
for a cut-off time-scale of several weeks; however, the uneven sampling of these data made their 
reliability uncertain. The observational situation was significantly improved with the launch of 
RXTE. Taking advantage of the unique properties of this satellite, Edelson and Nandra (1999) 
obtained evenly-sampled fluxes of the BLS1 NGC 3516 on long, medium and short time-scales, 
and combined these data to produce a PDS covering 4 decades in temporal frequency, finding a 
cut-off time-scale of ~1 month. 

Until now the only quantitative assessment of the rapid X-ray variability of NLSls has been 
based on measurements of the excess variance, obtained from ASCA observations, typically of 1 
day duration. Based on these data, Turner et al. (1999) and Leighly (1999) found the short-term 
variability of NLSls to be substantially greater than for BLSls, while in both cases the excess 
variance was anticorrelated with the X-ray luminosity. The present paper reports the results of the 
first extensive and evenly-sampled X-ray monitoring of an NLS1, using RXTE to monitor Akn 564 
over a near 2-year period. Akn 564 is particularly well suited for this study, being the brightest 
known NLS1 in the hard X-ray sky (F(2 — 10 keV) « 2 — 5 x 10~ n erg cm~ 2 sec -1 ; Vaughan et 
al. 1999a), strongly variable (e.g., 50% variations on timescales of hours; Leighly 1999) and located 
well out of the ecliptic plane (allowing monitoring by RXTE throughout the year). Subsequent 
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to the start of our RXTE monitoring campaign, Akn 564 was chosen for simultaneous monitoring 
with HST, ASCA, Chandra and the AGN Watch network of ground based optical telescopes. 

2. Observations and Data Reduction 

Our X-ray monitoring was planned to obtain variability information spanning time-scales from 
hours to months. Akn 564 was observed once every 4.3 days (= 64 orbits) from 1999 January 1 
- 2000 September 19 (the limit of the present analysis, although this low-frequency monitoring is 
continuing), and once every ~3.2 hr (= 2 orbits) from 2000 June 30 - August 1. (Twice-daily 
observations were also obtained during 1999 October 25 - November 11, in support of the AGN 
Watch campaign; Schemmer et al. in preparation). Our sampling parameters are summarized in 
Table 1. 

The RXTE Proportional Counter Array (PCA) consists of five collimated Proportional Counter 
Units (PCUs), sensitive in a nominal 2-60 keV bandpass (Jahoda et al. 1996). At the start of the 
present campaign three of the PCUs (0, 1 and 2) were in routine use. After the gain settings of 
all five PCUs on board RXTE were changed on 1999 March 22 this number was reduced to two 
(0 and 1), and then, from 2000 May 12, to one (1). Data collected during each of the two gain 
epochs were extracted separately using background model files appropriate to that gain epoch. 
The present analysis is restricted to the 2-10 keV band, where the PCA is most sensitive and the 
systematic errors are best understood. Data from the top (most sensitive) layer of the PCU array 
were extracted using the REX reduction script 4 , and background counts were estimated using the 
"L7-240" model 5 . Background-subtracted count rates were then obtained in each epoch to assign a 
mean flux and error. Count rates derived from data taken during the earlier gain setting were scaled 
to account for the effect of the gain change by a factor derived from observations of Cassiopeia A, 
and then normalised to units of ct s _1 PCU -1 . Figure 1 (top panel) reproduces the full light curve. 

3. Temporal Analysis 

Before performing the time series analyses described below, the observed light curve was re- 
sampled to provide two nearly-independent light curves, one sampling short time-scales (~days) 
and the other, long time-scales (~months). In Figure 1, the centre panel contains the full 20 months 
of data, sampled on a grid as close to 4.3 d as possible. Interpolating over the 9 data points that 
were missing during this period and mapping to an even grid yielded a total of 148 evenly-spaced 
flux points. This had the effect of filtering out (not smoothing) most of the data in the intensive 
periods around MJD 51480 and 51700. The short time-scale light curve, shown in the bottom 
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panel of Figure 1, was obtained during the 32 d intensive sampling period 2000 June 30 - August 1. 
Interpolating over the 13 data points missing during this period yielded a total of 235 evenly-spaced 
flux points. 



3.1. Comparison of Long and Short Timescale Light Curves 

The long and short time-scale data are quantitatively compared in Table 2, where column 2 
gives the mean count rates, and column 4 the fractional excess variance, cr^ s , defined as 

2 S 2 ~ i a err} /i\ 

°xs = {x)2 , (1) 

where (X) is the mean flux, (<J 2 rr ) is the mean square error, and S 2 is the measured variance of the 
light curve (e.g. Nandra et al. 1997). We have estimated (±1<t) errors on the excess variance as 
±ax S y/2/N, where N is the number of flux points, noting that for Gaussian fluctuations S 2 will 
follow a x 2 distribution (Warwick, private communication). 

The short time-scale light curve of Akn 564 shows a marginally larger excess variance than 
the long time-scale light curve, a result that is strikingly different from that of BLSls. In a survey 
of RXTE and ASCA data for 8 BLSls, Markowitz & Edelson (2001) found larger values of a 2 xs 
on long (~months) than on short (~hours) time-scales in every case. This comparison strongly 
suggests that large amplitude X-ray variability is more rapid, i.e., extends to significantly shorter 
time-scales, in Akn 564 than in BLSls of comparable luminosity. 



3.2. Fluctuation Power Density Spectrum 

In order to better quantify the short-term variability of Akn 564 the PDS was then constructed 
following the prescription of Edelson & Nandra (1999), first determining separate short and long 
time-scale PDS and then combining them to produce a single PDS. The individual PDS were 
derived using a direct Discrete Fourier Transform (Oppenheim & Shafer 1975, Brillinger 1981). 
The zero-power and next two (very noisy) lowest-frequency points of each PDS were ignored and 
the remaining points binned every factor of 1.7 (0.23 in the logarithm). Power-law models were then 
fitted separately to the long and short time-scale PDS. The slope and uncertainty, measured from 
an unweighted, least-squares fit to the logarithmically binned data, and derived using the Welch 
window function, are given in column 5 of Table 2. These show a highly significant (10<r) systematic 
flattening from short time-scales (with a = —0.96 ± 0.07) to long time-scales (a = —0.24 ± 0.08). 

The long and short time-scale PDS were then combined. Unlike Edelson & Nandra (1999), 
who allowed the relative normalization between PDS to float as a free parameter, the absolute 
normalizations were retained. The combined PDS is shown in Figure 2. A simple power law gave 
an unacceptable fit (>99.99%), with a reduced chi-squared of xl = 3-06 for 13 degrees of freedom. 



- 5 - 



The combined PDS was then fitted with a model in which a steep power law dominates at 
high frequencies, but cuts off to a slope of a = at low frequencies: 

P(/)=C/(1 + /// C ) a , (2) 

where P(f) is the fluctuation power at a frequency /, a is the power law slope at high temporal 
frequencies, f c is the "cutoff frequency," well below which the PDS flattens to a slope of zero, 
and C is the normalization. The fit, with three free parameters (C, f c , and a) was a significant 
improvement over the simple power law, with xt = 1-70 for 12 degrees of freedom, although 
still formally unacceptable (94% confidence) mainly due to the higher absolute level of variability 
power in the short time-scale data. The high-frequency slope of the broken power law fit is a = 
— 1.12, similar to the slope derived using the high-frequency PDS alone, and the cutoff frequency of 
8.7 x 10~ 7 Hz corresponds to a timescale of ~13 d. These parameters are summarized in Table 3. 
While the power law plus cutoff model is clearly a marked improvement over the unbroken power 
law, we note the exact shape of the turnover is not well constrained and the low-frequency behaviour 
not at all well determined. 



4. Discussion and Conclusions 

We have carried out the most comprehensive X-ray variability study of an NLS1 to date. 
The long and short time-scale light curves of Akn 564 are found to be remarkably similar, with 
comparable excess variances over the sampling intervals of 32 days and 627 days. The straight- 
forward deduction from a visual examination of the light curves is that Akn 564 exhibits most of 
its variability power on time-scales substantially less than one month. We suggest this may well 
be a common feature of NLSls, and one which distinguishes the sub-class from BLSls. Turner 
et al. (1999) calculated the excess variance for ASCA observations of 36 Seyfert galaxies over a 
typical sampling interval of 1 day finding values up to ~10 times greater for NLSls compared with 
BLSls (of similar luminosity). A current explanation for the rapid and large amplitude variability 
characteristic of Narrow Line AGN is in terms of a smaller black hole mass, and hence of size scale. 
Our result now confirms - for Akn 564 - that the large short-term variability is indeed primarily 
due to the time-scales being shorter (rather than - necessarily - having greater variability power) 
than for BLS1 of comparable luminosity. 

The Akn 564 PDS now offers an opportunity to quantify the lower mass and higher accretion 
rate of a NLS1 in comparison with a BLS1 of similar luminosity. To date the only BLS1 PDS 
obtained by a similar evenly-sampled campaign is for NGC 3516, from which Edelson & Nandra 
(1999) obtained a cut-off time-scale of ~30 d. Assuming for the moment a simple scaling law 
applies for the variability time-scale and black hole mass of NLS1 and BLS1, our result suggests 
that Akn 564 has a black hole mass of order ~ 0.4 that of NGC 3516. Reverberation mapping by 
Robinson et al. (1994) found NGC 3516 to have a black hole mass of order ~3x 10 7 M , implying 
a corresponding mass for Akn 564 of ~ 1 x 10 7 M . Estimation of the accretion rate for most 
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AGN is hampered by a large uncertainty of the luminosity in the hidden EUV band. However, 
Wandel et al. (1999) showed the bolometric luminosity of BLSls may be estimated by scaling from 
the continuum luminosity at 5100 A, based on a comparison of photoionisation and reverberation 
estimates of black hole masses. By this means we derive a bolometric luminosity for NGC 3516 
of ~ 10 44 erg/s which should be good to a factor of ~ 3. The corresponding accretion rate for 
NGC 3516 is then of order 0.05 M Ed d- 

The very different SED of NLSls means we cannot rely on the Wandel et al. (1999) relation to 
estimate the bolometric luminosity of Akn 564. Instead we use the respective 2-10keV luminosities, 
for which Akn 564 is some ~3 times that of NGC 3516. Combining the respective estimates of 
black hole mass and luminosity then gives an accretion rate some ~8 times larger for Akn 564 
than for NGC 3516, indicating that Akn 564 is accreting at a substantial fraction of the Eddington 
limit, or 0.4 MEdd- We note furthermore that this is probably a lower limit, given that the scaling 
with NGC 3516 via the 2-10 keV luminosities takes no account of the strong soft X-ray excess in 
Akn 564. 

Interpreting the PDS in the above way is subject to two caveats. First, is the variability 
stationary, to at least the extent that the measured " characteristic time-scale" is a reliable measure 
of scale and black hole mass? The uniformity of the long time-scale light curve of Akn 564 shown 
in Figure 1 is reassuring in that respect, suggesting the absence of low frequency power is sustained 
throughout our campaign. On the other hand, the higher mean flux and larger fractional variability 
during the intensive monitoring suggests a degree of non-stationarity, but not to an extent which 
alters our basic finding, that Akn 564 contains most of its variability power on timescales much 
shorter than 1 month. A second question is whether the accretion discs in both Akn 564 and 
NGC 3516 extend close to the innermost stable orbit, justifying our inverse scaling of black hole 
mass with the respective values of / c ? The high luminosity of NLSls indicate that is generally true 
for NLSls, while in the case of NGC 3516 the observation of a broad Fe-K fluorescence line (Nandra 
et al. 1999) suggests that its accretion disc also extends close to the innermost stable orbit. 

We conclude that comparisons of the X-ray variability of Akn 564 and NGC 3516 indicate a 
black hole mass for Akn 564 of ~ 1 x 10 7 M , implying an accretion rate in the range 0.2-1.0 
MEdd- It is interesting to compare this result with the accretion rate derived from modelling of line 
widths and variability time-scales in the optical Broad Line Region (Laor 1998). Taking the above 
estimate for the luminosity of Akn 564 and H/3 FWHM = 800 km/s (Vaughan et al. 1999b), the 
relation in Laor (1998) yields an "independent" estimate for Akn 564 of 1.0 MEdd 

The authors acknowledge the dedicated work of the RXTE team, especially Evan Smith for his 
careful scheduling so close to the ephemeris. RAE and AM were supported by NASA grants NAG 
5-7315 and NAG 5-9023, and SAV by the UK Particle Physics and Astronomy Research Council. 
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Table 1. Sampling Parameters 



Time 


Observing Dates 


Mean Sampling 


Data 


No. of Flux 


Usable Temporal 


Scale 


(MJD) 


Interval 


Lost 


Points 


Frequency Range 


Long 


51179.55-51806.78 


4.267 d 


9 points 


148 


6.9 x 1(T 8 - 1.1 x 1(T 6 


Short 


51694.82-51726.47 


3.2 hr 


13 points 


235 


1.4 x 10" 6 - 3.4 x 1(T 5 



Table 2. Multi-Band Light Curve Properties 



Time 


Mean 


Excess 


Power-Law 


Scale 


ct s- 1 


Variance 


Slope 


Long 


1.63 ± 0.04 


0.076 ± 0.009 


-0.24 ± 0.08 


Short 


1.87 ± 0.04 


0.111 ± 0.010 


-0.96 ± 0.07 



Table 3. PDS Fit Parameters 



Description 


Parameter 


value 


High-Frequency Slope 


a 


-1.12 ± 0.04 


Cutoff Frequency 


fc 


8.72 ±0.81 x 10" 7 Hz 


Cutoff Timescale 


tc(= l/fc) 


13.31^ days 


Normalization Coefficient 


Ci 


4.4 ±0.2 x 10 6 Hz" 1 
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Fig. 1. — Full band (2-10 keV) light curve (top), resampled on long (4.3 d-20 months; centre ) 
and short (3.2 hr-31 d; bottom) timescales. These middle and lower light curves were produced by 
filtering and interpolating the data in the top panel as described in the text. Interpolated points 
are shown as open circles without error bars, and the 8 points common to both light curves as x's. 
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Fig. 2. — The combined PDS of Akn 564. The crosses refer to the PDS derived from the short 
term light curve while the circles are from the long timescale PDS. The solid line is the best fit as 
described in the text, a power-law with a low-frequency cutoff at ~13 d. The dotted line is the 
same power-law without the cutoff. 



